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Calcium (Ca 2+ ) is a key secondary messenger in many plant signaling pathways. One such pathway is the SYM pathway, 
required in the establishment of both arbuscular mycorrhizal and rhizobial root symbioses with legume host plants. 1 
When the host plant has perceived the diffusible signals from the microbial symbionts, one of the earliest physiological 
responses are Ca 2+ oscillations in and around the nucleus. 2 These oscillations are essential for activating downstream 
gene expression, but the precise mechanisms of encoding and decoding the Ca 2+ signals are unclear and still under 
intense investigation. Here we put forward a hypothesis for the mechanism of the cation channel DM11. 



In a recent study 3 we presented a mathematical model based on 
three key components involved in the production of the Ca 2+ 
oscillations in the SYM pathway: a cation channel, 4,5 a Ca 2+ 
pump 2 and an as yet unidentified Ca 2+ channel. 1 The compo- 
nents involved in the model are located on the inner nuclear 
membrane 2 and therefore the model hypothesizes that the lumen 
of the nuclear envelope contiguous with the ER is the Ca 2+ 
store. The model also included Ca 2t -binding proteins, so called 
Ca 2 * buffers, which can be any Ca 2+ -binding proteins such as 
calmodulin, Ca 2 * /calmodulin dependent protein kinase 6,7 or 
Ca 2 * reporters used to measure Ca 2 * levels in biological systems. 8 
Our study revealed that varying the binding characteristics and 
concentrations of the buffers can affect both Ca 2+ spike shapes 
and period lengths, induce rapid Ca 2+ spiking and can initiate 
and terminate Ca 2 * oscillations. Overall, the mathematical model 
strikingly reproduces the Ca 2+ oscillations observed experimen- 
tally and suggests that buffering capacity, so far unexplored in 
this system, can explain changes in shape and frequency of the 
nuclear localized Ca 2 * oscillations in the SYM pathway. 

One of the three components used to generate the math- 
ematical model is the cation channel, named DMI1 in Medicago 
truncatula, which has been genetically identified and character- 
ized. 5,9 Although DMI1 is essential for generating nuclear Ca 2+ 
oscillations, its role is currently not well understood and prone to 
speculation. Upon activation by symbiotic secondary messengers, 



DMI1 could counter-balance the flow of positive charges gen- 
erated by the simultaneous activation of a yet to be unidenti- 
fied Ca 2 * channel. 4 Alternatively, the activation of DMI1 could 
directly trigger the opening of a voltage-gated Ca 2+ channel by 
hyperpolarizing the nuclear membrane potential. 5 These two 
hypotheses are experimentally unresolved, but we could use our 
mathematical model, that successfully recapitulates the experi- 
mental Ca 2 * oscillations, to test the potential functions of DMI1. 

As explained in our recent study, we envisage the membrane 
potential of the nuclear envelope to be negatively charged on the 
nucleoplasmic side relative to the perinuclear space. 3 In this con- 
text our mathematical model reveals that Ca 2+ spiking is only 
initiated upon the simultaneous activation of both DMI1 and 
the Ca 2 * channel; we hypothesize that this may be triggered by 
unidentified symbiotic secondary messengers. Once both chan- 
nels have opened, the generation and sustainability of the Ca 2 * 
oscillation is dependent on the interplay between the key play- 
ers DMI1 and the Ca 2+ channel via their respective K + and Ca 2+ 
electrochemical driving force and their conductance (Fig. 1). At 
the start of a spike (Fig. 1A and B) the membrane potential is 
negative and close to the K + resting potential, such that DMI1 
is weakly conducting (Fig. 1C, position 1). This negative mem- 
brane potential drives a Ca 2 * current into the nucleus (Fig. ID, 
position 1). The membrane potential depolarizes until it reaches 
the resting potential of Ca 2 * and consequently stops the transient 
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Figure 1. How DM11 controls the calcium oscillations. (A) An experimental trace of calcium oscillations determined by using microinjection of the 
dyes Oregon Green (OG) and Texas Red (TR) into a M. truncatula root hair cell. (B) A close-up of one of the spikes. Numbers (1,2,3,4) in Figures B-E are 
assigned to denote key stages or events during an oscillation, as explained in the main text. (C) Position 2-3, a calcium spike (red) begins at the same 
time as the membrane potential (black) decreases slowly. (D) Position 1 shows the first brief calcium current transient (red) that depolarizes the nu- 
clear membrane (as seen in the corresponding position 1, black, in Fig. C). Position 2-3, the calcium current (red) increases slowly as does the current 
through DM11 (black) in the opposite direction. (E) The figure illustrates the modulation of the calcium oscillations by DMIl.The color map shows how 
the potassium current released by DM11 depends on the calcium concentration and on the membrane potential; the larger the membrane potential is, 
the higher is the calcium concentration released and also the more distant the membrane potential is from the potassium resting potential (-17.7 mV). 
Superimposed on the color map, in black dots, we show the phase space diagram of Ca 2+ concentration and voltage. The system trajectory oscillates 
along with the potassium current: the calcium concentration rises when the magnitude of the potassium current is increasing (branch 2), and it falls 
when the potassium current is very low (branch 4, and corresponding Position 4 in Fig. B-D). The parameters are the same as in Granqvist et al. 3 



Ca 2t current (Fig. 1C and D, position 1). This first step has two 
important consequences: the localized release of Ca 2+ increases 
the opening of DMI1, and the positive membrane potential 
increases the K + electrochemical driving force. 

Subsequently the current generated by DMI1 slowly increases 
and hyperpolarizes the membrane potential, as K + flows into the 



perinuclear space (Fig. 1C and D, position 2). This hyperpolar- 
ization generates a Ca 2+ current (Fig. ID, position 2-3) leading 
to the Ca 2+ release which shapes the upward slope of a Ca 2+ spike 
(Fig. 1C and D, position 2-3). As soon as the membrane poten- 
tial reaches the resting potential of K*, the DMI1 current and 
the Ca 2+ channel current almost cease (Fig. ID, position 3-4). 
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The Ca 2 * is pumped back into the store by Ca 2 *ATPases such as 
MCA8, 2 decreasing the nucleoplasmic Ca 2+ concentration which 
shapes the downward slope of a Ca 2 * spike (Fig. IB, position 
4). Although both channels are weakly conducting, the electro- 
chemical driving force of the Ca 2+ channel is very strong (Fig. 
IE, position 4). As soon as the membrane potential returns to 
the initial value, the conductance of the depolarization-activated 
Ca 2 * channel increases again leading to a release of Ca 2 *. The 
Ca 2 * release increases DMI1 conductance and subsequently the 
DMI1 current hyperpolarizes the membrane; the cycle starts and 
repeats to generate a Ca 2+ oscillation. 

Overall, the mathematical model suggests that DMI1 acts 
predominantly as a counter ion channel, but importantly initial 
conductance from DMI1 is necessary to facilitate calcium flow. 
In this regard DMI1 function is neither purely a counter ion 
balance, nor purely an activator of the calcium channel. Rather 
the interplay of DMI1 and the calcium channel together derive 
the oscillatory calcium behavior. The model predicts that a key 
parameter of DMI1 action is modulation by Ca 2 * in order to 
sustain an oscillatory mechanism. Indeed the increased K + cur- 
rent generated by DMI1 synchronously triggers the Ca 2 * release 
(Fig. IE) and can be obtained by direct positive feedback of Ca 2+ 
on DMI1 conductance. This positive feedback could potentially 
occur by direct binding of Ca 2 * to DMI1. In agreement with 
this hypothesis, a previous study suggested the presence of Ca 2 *- 
binding pockets in the C-terminal region of DMI1 based on 
homology modeling with its closest structural relative MthK. 10 In 
addition, Venkateshwaran and colleagues 5 recently demonstrated 
that expression of DMI1 in HEK cells is sufficient to activate 
Ca 2 * induced Ca 2 * release mechanism upon Ca 2+ stimulation. 
Therefore a positive Ca 2+ feedback could explain the oscillatory 
mechanism but also the sustainability of the Ca 2 * oscillation. 
Indeed, the oscillation is sustained until the buffering capacity of 
the cell breaks the positive Ca 2+ feedback. Ca 2+ could be buffered 
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such that the Ca 2+ concentration is no longer sufficient to modu- 
late DMI1 conductance and would therefore terminate the oscil- 
lation. 3 The positive feedback mechanism of Ca 2 * on DMI1 
conductance however has yet to be experimentally validated. 

A recent study suggested that the role of DMI1 is to activate 
a yet unidentified voltage-gated Ca 2+ channel by hyperpolariz- 
ing the nuclear membrane potential, and that its ability to do so 
depends on the length of DMI1 open time. 5 We challenged our 
model with these hypotheses and could only generate one Ca 2 * 
spike but no Ca 2 * oscillations. In this hypothetical scenario, a 
sustainable Ca 2+ oscillation would only be obtained if the con- 
centration of the secondary messenger activating DMI1 would 
oscillate synchronously with the Ca 2 * transient. This speculation 
brings an additional level of complexity to the generation of the 
Ca 2 * oscillation, which has never been reported so far in any Ca 2 * 
induced Ca 2 * release mechanism. 

In summary, our model suggests that DMI1 plays a key role 
as a counter ion channel in the generation of Ca 2 * oscillations. 
However, this counter ion channel function is intrinsic to the 
generation of the calcium current and initial K* movement facili- 
tates the first steps of calcium release. Without the K* current no 
Ca 2 * oscillation can be sustained. To validate this mathematical 
model, future experimental work is required, notably to identify 
the symbiotic secondary messengers activating the channels, 
define the effect of calcium on DMI1 function and to character- 
ize the missing player, the Ca 2 * channel. 
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